Drude weight, cyclotron resonance, and the Dicke model of graphene cavity QED 
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The Dicke model of cavity quantum electrodynamics is approximately realized in condensed matter 
when the cyclotron transition of a two-dimensional electron gas is nearly resonant with a cavity 
photon mode. We point out that in the strong coupling limit the Dicke model of cavity cyclotron 
resonance must be supplemented by a term that is quadratic in the cavity photon field and suppresses 
the model's transition to a super-radiant state. We develop the theory of graphene cavity cyclotron 
resonance and show that the quadratic term, which is absent in graphene's low-energy Dirac model 
Hamiltonian, is in this case dynamically generated by virtual inter-band transitions. 



Introduction — Distinct but degenerate dipole transitions 
between matter states can behave cooperatively in an 
optical cavity because the photon continuum is dis- 
cretized pQ. This physics is often captured approximately 
by the Dicke model [2] which describes a non-dissipative 
closed system of identical two-level systems interacting 
with a single-mode radiation field. For a sufficiently 
strong light-matter coupling constant, the thermody- 
namic limit of the Dicke model exhibits a second-order 
quantum phase transition to a super-radiant ground state 
(SPT) [3 with macroscopic photon occupation and co- 
herent atomic polarization. Wide interest in these SPTs 
has emerged recently in the context of circuit quantum 
electrodynamics (QED) 4J. It has been shown in particu- 
lar that superconducting flux qubits inductively coupled 
to a transmission line resonator display a SPT charac- 
terized by a doubly-degenerate ground state with macro- 
scopic photon occupation [5]. There is also interest in 
the possibility of a related phase transition for Cooper- 
pair boxes capacitively coupled to a transmission line 
resonator [51 [7]- Ultracold atom gases in optical cavi- 
ties [H [pj , in which spontaneous breaking of translational 
symmetry has been reported [5], have provided still an- 
other interesting new application of Dicke model physics. 

When an external magnetic field is applied to a two- 
dimensional (2D) electron system, transitions between 
states in full and empty Landau levels (LLs) are dis- 
persionless, mimicking atomic transitions and enabling a 
condensed matter realization of the Dicke model. In par- 
ticular, recent pioneering theoretical [10] and experimen- 
tal 1 11 j work has shown that these systems can be driven 
toward the ultrastrong coupling 12 limit by tuning the 
cyclotron transition energy of an ordinary parabolic-band 
2D electron gas to resonance with the photonic modes of 
a terahertz metamaterial. 

The light-matter interaction in the Dicke Hamiltonian 
is linear in the vector potential A cm of the cavity. For 
condensed matter states described by parabolic band 
models, a quadratic A^. m term whose strength is related 
to the system's Drude weight [13] and f-sum rule [14] , also 
emerges naturally from minimal coupling. It has long 
been understood jT5] that the Dicke model's SPT is sup- 



pressed when the quadratic terms are retained. Demon- 
strations of this property are often referred to as no-go 
theorems. (Standard no-go theorems do not apply |16j 
to ultracold atoms which are driven by an external pump 
field and subject to significant cavity losses.) 

We focus here on the case of cyclotron transitions in 
a graphene 2D electron system. Electronic states near 
the neutrality point of a graphene sheet [T7] exhibit lin- 
ear crossing between conduction and valence bands. The 
effective low-energy theory for electrons in graphene has 
a 2D massless Dirac fermion (MDF) Hamiltonian [18) . 
which is linear in momentum p. As first noted in 
Ref. [IH] , minimal coupling applied to the MDF Hamilto- 
nian does not generate a term proportional to A^ m . Cy- 
clotron resonance in this material, which has been exten- 
sively investigated experimentally and theoretically over 
the past decade [2UJ , seems therefore to provide an exam- 
ple of an active medium which could enable a SPT [TpJ 
when the graphene sheet is embedded in a cavity. In- 
deed, recent experimental advances have made it pos- 
sible to monolithically integrate graphene with optical 
microcavities [2TJ [22] , paving the way for cavity QED at 
the nanometer scale with graphene as an active medium. 
In this work we demonstrate however that in the strong 
coupling regime the Dicke model for graphene cavity cy- 
clotron resonance must be supplemented by a quadratic 
term that is dynamically generated by inter-band tran- 
sitions and again implies a no-go theorem. These no-go 
theorems are ultimately a consequence of gauge invari- 
ance and apply whenever this symmetry is unbroken. 
Gauge invariance and SPTs — We consider an electronic 
system in D spatial dimensions coupled to an electro- 
magnetic (e.m.) field with a single privileged mode de- 
scribed by a vector potential A(r,t). We argue below 
that the no-go theorem for the SPT requires only unbro- 
ken gauge symmetry. Our ideas are most clearly spelled 
out when A is treated classically. Quantization of the 
e.m. field can be easily carried out in the final step. 
Light-matter interactions are described by minimal cou- 
pling: pi — > pi + eA(ri,t)/c, where Pi is the canonical 
momentum of the z-th electron and — e is the electron 
charge. The Hamiltonian of a light-matter system can 
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always be written as H[A] = V. ma t[A] + % cm [A], where 
H ma t[A] contains all the electronic degrees of freedom 
treated quantum mechanically, while K cm [A] is the clas- 
sical energy density of the cm. field. The spontaneous 
coherent photon state is the ground state when the to- 
tal energy is lowered by introducing a finite static vector 
potential. Since 'Hem [A] is a positive-definite quadratic 
form of A, the instability can occur only if the second 
derivative of the matter energy with respect to A is neg- 
ative for static A. 

We therefore consider the variation of the matter en- 
ergy AE mat = E mat [SA]-E mat [0] due to an infinitesimal 
variation of the static vector potential: 



undoped 



doped 



J d D r 5A-(SH mat [A]/SA) 



(1) 



The quantity 5'H ma t[A}/SA is [33] the physical current 
operator, «7 p hy S (^)- It is convenient in perturbative anal- 
yses to split the current into paramagnetic and diamag- 
netic contributions, which can be defined by the expan- 
sion of J p hys( r ) m powers of 5 A. For the /z-th Cartesian 
component of the current operator this expansion reads 



J$ t) (r) + ~ J d D r'K^(r : r')SA„{r') 



(2) 



where J^\r) = <F# ma t[A]/&A At (r)| a=o is the para- 
magnetic current-density operator, K ,iv {r,r') — 
(c/e)S 2 H mat [A]/8A ll (r)6A u (r')\ A=0 , and the sum over 
repeated Greek indices is intended. The second term on 
the r.h.s. of Eq. ([2| is the diamagnetic contribution to 
the current-density operator. 

We now evaluate the expectation value of Jp^ ys (r) to 
first order in 8 A: i) the expectation value of the para- 
magnetic current-density operator, Jg(r), can be cal- 
culated using linear response theory [131 I14j : we find 
J£(r) = V~ 1 J2 q Jp(<l)e iq ' r + cc. with V = L D the 
/^-dimensional volume and Jg(q) = (e/ 'c)x Mt '(q)SA u (q), 
where the tensor x^il) is the static paramagnetic 
current-current response function. This current corre- 
sponds to the linear term in energy that is retained 
in the Dicke model. The diamagnetic current operator 
in Eq. ^ corresponds to the quadratic term neglected 
in the Dicke model. It already contains one power of 
5 A. Its contribution to the linear response current can 
therefore be obtained by evaluating its expectation value, 
K^(\r~r'\) ee (k» u {r,r')), where (...) refers now to 
the ground state of H ma t at 5 A = [24] . 

Using these results, we find that (t/phy S ( r )) = 
(e/cjV^ 1 ^2 q E^" (q)5A„(q)e lq " r +c.c, where we have in- 
troduced E^(q) ee x ^{q) + K^(q), K^(q) being the 
Fourier transform of K^ v {\r — r'\). The corresponding 
change in matter energy up to second order in 8 A is 
AE mat cx V- 1 J2^(q)SA fl (q)5A l/ (~q). The key point 
now is to realize [131 114) that, if gauge invariance is un- 
broken, there cannot be a current response to a static 
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FIG. 1: (color online) Left panel: Dipole allowed transitions 
in a neutral graphene sheet in the presence of an external 
magnetic field. The horizontal lines denote the unevenly- 
spaced Landau levels of massless Dirac fermions. The filled 
valence band levels (band label A = — 1) and the empty con- 
duction band levels (band label A = +1) are indicated by 
filled (empty) circles. The zero-energy (n — 0) Landau level 
is formed partly from the valence band and partly from the 
conduction band and is half filled in a neutral system. Right 
panel: In a doped graphene sheet with a Fermi level ef placed 
between the conduction band n — M and the n = M + 1 
Landau levels, there is an allowed cyclotron transition within 
the conduction band (magenta line) while some of the lower 
energy inter-band transitions, indicated by dashed lines, are 
Pauli-blocked. There is a clear energetic separation between 
the lowest energy (unblocked) intra-band transition and the 
(unblocked) inter-band transitions. 



uniform potential [25) . The diamagnetic contribution to 
the current- current response is crucial to cancel the para- 
magnetic contribution and can never be neglected: 



lim K^(q) 

q-s-0 



lim X ^(g) 

q-*0 



(3) 



This result, which is the most important of this Section, 
implies that, in the absence of broken gauge invariance, 
the change in matter energy due to a static and quasi- 
homogeneous vector potential must vanish |26j . SPTs 
therefore cannot occur unless gauge invariance is broken, 
as occurs for example in the case of superconductors. 
Cavity cyclotron resonance in graphene — We now focus 
on a 2D electron gas in a graphene sheet exposed to a 
strong quantizing magnetic field B = Bz [29]. When 
Zeeman coupling is neglected, the LLs of graphene have 
two-fold spin and valley degeneracies and therefore over- 
all degeneracy Nf — 4 [5D]. The single-spin single- 
valley MDF Hamiltonian reads Ho = vr>cr ■ II(r), where 
vr> fa 10 6 m/s is the Dirac velocity, er = (a x ,a v ) is a 
2D vector of Pauli matrices, H(r) = —ihV r + eA (r)/c 
is the gauge-invariant dynamical momentum, and A (r) 
is the static vector potential that describes the quan- 
tizing magnetic field. We work in the Landau gauge 
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A (r) = (—By, 0,0). A complete set of eigenfunctions 
is provided by the two-component pseudospinors 20 



{r\X,n,k) 



L 



c n (p„-i,k(y) 

\C+(p n>k (y) 



(4) 

where A = +1 ( — 1) for conduction (valence) band levels, 
n E N is the intra-band LL index, and k is the eigen- 
value of the magnetic translation operator in the x direc- 



tion. In Eq. (|4j) = y^l ± 5 n ^ )/2 and <f n ,k(y) are the 
normalized eigenfunctions of the parabolic-band Landau 
problem [TJ] . The form of the coefficients guarantees 
that the pseudospinor corresponding to the n = LL has 
weight only on one sublattice. Each LL has a degener- 
acy Af — N{S/ (2tt£ 2 b ), where S — 1? is the sample area. 
The pseudospinor ([4]) is an eigenstate of the Hamiltonian 



-(0) _ 



= \huj r 



where we have 



Ho with eigenvalue [^y c A , 
introduced the MDF cyclotron frequency lu c — vd\/2/£b 
and the magnetic length is = [he/ {eB)] 1 / 2 . 

The Hamiltonian which describes coupling between 
MDFs and light in the cavity does not contain a quadratic 
term: Hint = vu(e/c)cr ■ A em (r). For future purposes 
we introduce the notations a ± — (a x ± ia v )/2 and 
^em( r ) = A^ m (r) ± zA^ m (r). In what follows we ne- 
glect [115] the spatial variation of the e.m. field in the 
cavity, Af m (r) — > A^ nl , since the photon wavelength is 
normally much larger than other length scales in the 
problem. In this quasi-uniform approximation we can 
easily evaluate the matrix elements of Hint between the 
unperturbed pseudospinors Q: 



(A ,n ,fc |W in t|A,n,fe) 



evr> 



n' ,ro-|-l 



+ A'C+,C n 5 n ' )n _iA+ m ) 



(5) 



The strong coupling limit is most easily obtained when 
the Fermi energy e-p lies within one of the bands; we con- 
sider the case in which it lies in the conduction band (A = 
+1) between the LL with index n = M, which is fully oc- 
cupied, and the LL with index n = M+l, which is at least 
partially empty. (See Fig. [I]) The Dicke model of cavity 
cyclotron resonance includes only the intra-band n = M 
to n = M + 1 transition and acts in the 2-^-fold subspace 
spanned by {\+,M,k), \+,M + l,k),k = 1...AT}, ne- 
glecting inter-band transitions [19] . Using Eq. ^ and in- 
troducing a set of Pauli matrices {1^, r|, r fe , k = 1 . . .Af} 
which act in this two-level-system subspace leads to the 
following pseudospin Hamiltonian: 



E 



n 



M 



~ T k 



c cm ' k 



(6) 



where Em 



1 - y/M 



M), and e e 



= evY>A~ m /(2c). In this 



model the occupied conduction-band LL shifts down in 



energy by 



AE 



(intra) 
M 



2 At 



l -(VM + 1 + VM) (7) 



in the limit of a static vector potential, in violation of 
gauge invariance as explained in the previous section. 
This is the origin of the SPT found in Ref. [Hj. The 
correct effective matter Hamiltonian 7i mat for graphene 
cavity cyclotron resonance must repair this defect. 

The Dicke model misses a diamagnetic contribution to 
7i mat , which, according to (J3j , must precisely cancel the 
spurious energy shift ((7j). To derive this term wc first 
recognize the intrinsic two-band nature of graphene (see 
Fig. [I]). A generic valence band state |— ,n, k) is cou- 
pled by the e.m. field to two states in conduction band: 
|+,n + 1, k) and \+,n — 1, k). We first consider the un- 
doped limit in which all valence band states |— , n, k) are 
occupied. Because the Dirac model applies over a large 
but finite energy region we must apply a cut-off by occu- 
pying valence band levels with < n < v max . Treating 
the e.m. field again by second-order perturbation the- 
ory, we find the following change in matter energy for an 
undoped graphene sheet in a quantizing magnetic field: 



AE 



undoped 



Af i/ max 

EE 



k=l n=0 



\(+,n+l,k\H iu t\-,n,k)\' 
JO) jo) 

t -,n t + .n+l 



Pn-1- 



|(+,n - l,fc|H int |- n,fc)| ; 

JO) _ JO) 



, (8) 



where p n — 1 — 8 n ^/2. (The factor p n takes care of 
transitions involving the n = LL, which is half filled.) 
Using Eq. ^ for the matrix elements we can write Eq. ^ 
more explicitly: 



AE, 



undoped 



(9) 



with F(u) = ELi (; 



- , — r - r i - Once 
again, this large negative contribution to the change in 
matter energy is spurious. It is present because the Dirac 
model with a rigid ultraviolet cut-off v max breaks gauge 
invariance [57] ■ When a model that is correct at atomic 
length scales, for example a 7r-band tight-binding model, 
is used instead, a static vector potential merely reassigns 
momentum labels within the full valence band. We com- 
pensate exactly for this deficiency of the Dirac model at 
its ultraviolet cut-off scale by adding the positive quan- 
tity — Ai? un( ioped to the change in matter energy. 

We now reconsider the situation analyzed earlier in 
which the Fermi energy ep lies in conduction band (A = 
+1) between LLs with indices n = M and n = M + 1, 
but account for inter-band transitions. The inter-band 
correction to the energy shift Ai?j^ ntra ' ) in Eq. (H), can 
be calculated using an expression which is equivalent to 
Eq. ^ but accounts for Pauli blocking of transitions to 
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occupied conduction-band states. The final result for the 
inter-band contribution is given by 



AE 



(inter) 
M 



-AE, 



undoped 



2 A 2 



hbJ c 



x - F{M)} , 

= N (^) 2 1k [1 + FiM)] > (10) 



where the term — A£? unc ;opcd takes care of the Dirac 
model regularization and Eq. ^ has been used in the last 
equality. After noticing that 1 + F(M) = sjM + 1 + a/M, 
we finally obtain 

A< ntcr) =Af(^y ^k(VM + l + VM) . (11) 

The quantity AE^ is a dynamically-generated inter- 
band diamagnetic contribution to the effective Hamilto- 
nian T-L e s, which: i) is independent of cut-off u max and ii) 



satisfies AE ma , t — AE^ tTa ^ 



A£; (intcr) = ^ . e . n 



limit of a static vector potential, it precisely cancels the 
spurious shift ^ responsible for the Dicke model SPT. 

Because intra-band transition energies are much lower 
than inter-band transitions in the weak-field limit rel- 
evant to the strong coupling limit of cavity cyclotron 
resonance, we can neglect the frequency dependence of 
the dynamically generated quadratic term. This energy 
must be added to the effective matter Hamiltonian Q 
for graphene cavity cyclotron resonance: 



H 



eS 



H 



eS 



s Vm A 2 



(12) 



where T>m = 4£M&uni/fi is the Drude weight ex- 
pressed in terms of the function £m introduced right af- 
ter Eq. ^ . This Hamiltonian is the starting point of the 
cavity QED theory of graphene cyclotron resonance. The 
A 2 m quadratic supplement to the Dicke model is always 
critical in the strong coupling limit. Eq. (12 1 is the most 
important result of this work. 

Quantum Theory — We can quantize the e.m. field by pro- 
moting the positive and negative Fourier amplitudes of 
A cm to photon annihilation a and creation operators: 
^4-cm = Ae(a + aT), where e is unit vector describing the 
polarization of the e.m. field and A = he 2 / (euiV) 
with V — SL Z the volume of the cavity (L z <C L is 
the height of the cavity in the direction perpendicular to 
graphene) and e its dielectric constant. When a cavity 
model with frequency u> is nearly resonant with the cy- 
clotron transition frequency Qm, the total Hamiltonian 



( 12 ) yields a Dicke model supplemented by an A 2 m term: 



u 



Dicke 



fyjja* a 



fc=i 



Af 



k=l 



+ K(a + a f ) 2 



(13) 



where g = hui c 2cr uni /(eu>L z ) and k — KDj^/(eu>L z ). In 



writing Eq. ( 13 1 we have assumed a specific polarization 
of the e.m. field, i.e. e = x. In the thermodynamic 
J\f — » oo limit the model (131 undergoes a SPT if the 
condition wO M (l + 4«;/w)/(4g 2 ) < 1 is satisfied [71 [T5]. 
In our case, however, a SPT is forbidden because the 
following identity holds true: 
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M 



(14) 



Equation ( 14) specifically establishes a no-go theorem for 
the occurrence of a SPT in the graphene cyclotron res- 
onance cavity QED. It is not a coincidence, and instead 
follows directly from the cancellation between paramag- 
netic and diamagnetic currents discussed in the first part 
of this work. The paramagnetic response of the Hamilto- 
nian ( |13[ ) to a static and quasi-homogeneous e.m. field is, 
indeed, (g 2 /N) lim^o((7f ot ; rf ot )) w = -2g 2 /n M , where 
T-tot = E k rt and (ferfj). = 2Afn M /(u 2 -tf M ) 



According to Eq. ([3]), this paramagnetic contribution 
must be equal in magnitude and opposite in sign to the 
diamagnetic response of (13), which is simply 2k: i.e. it 



must satisfy k — g /Qm, which coincides with Eq. (14). 
Summary — We have derived a microscopic effective 
Hamiltonian - Eq. ( 12 ) - for graphene cavity cyclotron 
resonance, highlighting in particular the role of gauge 
invariance and Drude weight. We find that the model 
must include a quadratic term which eliminates the pos- 
sibility of super-radiant quantum phase transitions. The 
peculiar spectrum of the massless Dirac fermion Hamil- 
tonian in a quantizing magnetic field, in which the Lan- 
dau level spacing decreases with increasing energy, nev- 
ertheless makes this material particularly attractive for 
ultrastrong coupling between cavity photon modes and 
cyclotron transitions. Electrons in graphene are thus 
an extremely intriguing active medium, and might pave 
the way for the realization of a gate-tunable general- 
ized Dicke model. Light-matter interactions can be fur- 
ther enhanced by combining cavity-integrated graphene 
sheets with graphene-based plasmonic elements |31j . Our 
work does not touch upon the interesting possibility of 
realizing super-radiant quantum phases when massless 
Dirac fermions are driven far away from equilibrium. 
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